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EXECUTIVE  SUMMARY 


This  final  report  presents  research  carried  out  by  Berkeley 
Research  Associates,  Inc.  under  Contract  F4 9620-8 6-C-0076  with 
the  Air  Force  Office  of  Scientific  Research  (AFOSR) .  The  report 
covers  the  performance  period  June  1986  through  February  1988. 
During  this  period  we  carried  out  research  on  the  modeling  of 
atomic  processes  for  x-ray  laser  plasmas.  The  report  is  divided 
into  two  parts.  Part  I  covers  research  performed  from  June  1986 
through  July  1987  and  Part  II  covers  research  performed  from 
August  1987  through  Februar  1988. 

In  Part  I,  we  deal  with  the  theoretical  modelincr  and 
computation  of  photoionization,  photoexcitation,  radiative  decay 
and  radiative  recombination  processes  for  multielectron  atoms  and 
ions  in  plasmas.  Such  atomic  processes  lead  to  and  influence 
lasing  in  x-ray  lasers  operating  both  in  the  soft  and  hard  x-ray 
regions.  Most  previous  calculations  for  multielectron  ions  have 
used  a  single  electron  or  independent  particle  model.  Although 
comparison  of  experimental  data  with  such  models  shows  agreement 
for  simple  systems  such  as  a  neutral  few-electron  atom, 
substantial  discrepancies  are  found  for  many  electron  atoms  and 
ions  with  a  large  number  of  bound  electrons.  In  our  research  we 
utilized  a  two-component,  finite  temperature,  self  consistent 
density  functional  method  and  demonstrate  that  this  method  is 
applicable  for  arbitrary  plasma  density  and  temperature  and  is 
capable  of  accurately  treating  multielectron  ions  of  arbitrary  Z. 

In  Part  II,  we  investigate  electron  collisional  ionization 
and  excitation  processes.  These  processes  are  an  important 
mechanism  through  which  population  inversion  of  ionic  energy 
levels  lead  to  lasing  in  the  soft  x-ray  region.  Previously 
theoretical  calculations  utilizing  a  semi-classical  impact 
approximation  have  been  performed.  In  many  cases  such  a  method 
is  inadequate.  In  our  research,  we  carried  out  computation  of 
electron  impact  ionization  and  excitation  cross-sections  and  rate 
coefficients  utilizing  the  distorted  wave  with  exchange  method. 

We  found  that  this  approach  provides  a  speedy  and  reasonably 
accurate  tool  for  these  calculations. 


PART  I 

MODELING  OF  ATOMIC  PROCESSES 
FOR 

X-RAY  LASER  PLASMAS 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  1 

II.  THE  METHOD  OF  CALCULATION . 6 

a.  Isolated  Atoms  or  Multielectron  lens  .  6 

b.  Partial  Cross-Sections  .  9 

c.  Radiative  Decay . 11 

d.  Radiative  Recombination . 11 

e.  Inclusion  of  the  Effect  of  Plasma 

Density  and  Temperature . 11 

III.  RESULTS . 16 

IV.  SUMMARY . 22 

V.  ACKNOWLEDGEMENTS  .  25 

VI.  REFERENCES . 2  6 

VII.  TABLES . 27 

VIII.  FIGURES . 31 


I .  INTRODUCTION 

Development  of  X-ray  lasers  operating  both  in  the  soft 
and  hard  X-ray  region  is  currently  an  area  of  intensive 
research.  A  variety  of  atomic  processes  lead  to  and 
influence  lasing  action  in  such  systems.  Photoionization, 
photoexcitation,  radiative  decay  and  radiative  recombination 
are  some  of  these  processes.  This  report  deals  with  the 
theoretical  modeling  and  computation  of  the  above  mentioned 
processes  for  multielectron  atoms  and  ions  in  plasmas. 

Accurate  calculation  of  photoionization  and 
photoexcitation  cross-sections  of  atoms  and  ions  are  useful 
in  a  variety  of  investigations  in  plasma  physics  and  atomic 
physics.  It  is  particularly  useful  in  the  context  of 
photopumping  schemes  for  X-ray  lasers.  As  another  example, 
computation  of  opacities  of  plasmas  for  diagnostic  and 
target  response  effects  require  these  data  as  input.  In 
order  to  model  the  radiation  spectra  from  hot  plasmas  (via 
the  detailed  configuration  rate  equation  for  example),  these 
data  are  required  along  with  the  other  bound-bound,  bound- 
free  and  free-free  processes.  Radiative  recombination  and 
radiative  decay  rates  are  very  useful  in  X-ray  laser 
research  -  particularly  in  the  context  of  recombination 
schemes .  In  view  of  these  different  applications,  there  is 
a  need  for  realistic  modeling  of  these  processes  to  generate 
accurate  data  for  a  variety  of  atoms  and  ions. 
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In  a  large  number  of  existing  calculations,  the  above 
processes  have  been  calculated  using  simple  approximate 
formulae,  valid  for  hydrogenic  ions  only.  Depending  on 
plasma  conditions,  however,  multielectron  ions  (Ne-like  ion, 
for  example,  with  ten  electrons  bound  to  it)  may  have  a  much 
higher  abundance  than  hydrogenic  ions.  Most  calculations 
for  multielectron  ions  existing  in  the  literature  have  been 
done  using  the  single  electron  or  independent  particle  model 
(IPA).  In  this  model,  the  energy  levels  and  wavefunctions 
of  the  atom  or  ion  are  first  calculated  using  the  Hartree- 
Fock  (HF)  method.  For  photoionization  calculation,  the 
interaction  of  the  incident  electromagnetic  radiation  with 
the  atom  (or  ion)  is  treated  via  the  first  order 
perturbation  theory1. 

Comparison  of  experimental  data  with  the  IPA 
calculations  shows  that  for  some  simple  systems  such  as  a 
neutral  few-electron  atom  (Lithium,  for  example),  there  is 
qualitative,  sometimes  quantitative  agreement.  However,  for 
many  electron  atoms  (Xenon,  for  example)  and  ions  with  a 
large  number  of  bound  electrons,  substantial  discrepencies 
are  found  between  experimental  and  iPA-data1. 

In  our  work,  we  used  the  time  dependent  linear 
response  approximation  within  the  framework  of  relativistic 
local  density  functional  method  ( DFM ) ^ '  ’  to  treat  the 
problem  of  photoionization.  This  method  incorporates 
several  advantages  over  the  HF-method.  The  HF-method  is 
non-local  and  computationally  very  elaborate,  whereas  in  the 
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density  functional  method,  one  deals  with  a  set  of  local 

equations  only.  This  leads  to  computational  simplicity.  On 

the  other  hand,  from  extensive  applications  of  the  density 

functional  method  (DFM)  it  is  well  known  that  fairly 

accurate  atomic  energy  levels,  wavefunctions  etc.  are 
1 4 

obtained  .  The  computational  simplicity  is  even  more 
apparent  in  the  case  of  relativistic  DFM  versus  relativistic 
HF-methods,  (proper  treatment  of  high  Z  atoms  require  the 
application  of  the  relativistic  DFM).  In  the  DFM, 
correlation  effects  of  the  bound  electrons  in  the  atom  are 
accounted  for  in  a  simple  way  via  the  correlation 
potential^.  The  HF-method,  on  the  other  hand,  does  not  take 
into  account  electron  correlation,  although  it  accounts  for 
non-local  exchange  effects  appropriately. 

For  photoionization  and  photoexcitation  process  (and 
therefore  for  radiative  recombination,  since  this  is 
effectively  the  inverse  of  photoionization),  the 
polarization  effect  of  the  atom  brought  about  by  the 
incident  radiation  field  is  important.  The  independent 
particle  method  does  not  take  this  into  account.  In  the 
linear  response  method  within  the  DFM  we  have  used  in  this 
work,  this  effect  is  treated  adequately  -  as  will  be  seen 
from  comparison  with  experimental  data.  In  most 
experimental  situations,  the  incident  radiation  (from 
synchrotron  sources  or  lasers)  have  field  strengths  small 
compared  to  the  atomic  field  strengths.  For  those 
experimental  conditions,  the  present  model  based  on  linear 
response  is  adequate  and  useful. 
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The  above  discussion  deals  with  the  calculation  of 
atomic  structure  and  processes  for  an  isolated  atom  or  ion  - 
that  is,  without  the  effect  of  the  plasma  environment.  This 
is  valid  for  very  low  density  plasmas.  As  the  plasma 
density  increases,  several  additional  effects  such  as  the 
screening  shifts  of  energy  levels  (or  lowering  of  ionization 
potential)  and  modification  of  wave-functions  and  potentials 
for  the  ion  embedded  in  high  density  plasmas  have  to  be 
considered.  The  effect  of  plasma  density  on  the  atomic  (or 
ionic)  structure  and  atomic  processes  needs  to  be  properly 
investigated.  In  the  present  work,  we  investigated  these 
effects  via  the  two-component  finite  temperature  self- 
consistent  DFM^  and  compared  the  results  with  corresponding 
isolated  ion  calculations. 

The  study  of  plasma  density  (and  temperature)  effect 
is  also  important  to  determine  the  optimum  operating 
conditions  for  gain  enhancement  in  X-ray  lasers.  At  very 
low  plasma  density,  collisional  excitation  and  deexcitation 
rates  are  small;  radiative  decay  is  then  a  dominant  process. 
However,  as  the  plasma  density  increases,  the  collisional 
rates  increase  rapidly  and  collisional  radiative  equilibrium 
can  exist.  These  in  turn,  have  their  effects  on  the 
relative  abundance  of  ionic  population  as  well  as  on  the 
population  inversion  of  ionic  levels  and  thus  affect  lasing. 
The  investigation  of  plasma  density  and  temperature  effects 
on  atomic  processes  is  therefore,  of  importance  in  overall 
X-ray  laser  modeling. 
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The  self-consistent  two  component  DFM  method  we  have 
used  provides  a  useful  tool  for  realistic  modeling  of  the 
various  effects  arising  in  high  density  hot  plasmas.  The 
electron  degeneracy  effect  is  correctly  treated  via  the 
Fermi  distribution  function  and  the  exchange-correlation 
effect  of  bound  and  free  electrons  are  taken  into  account 
via  the  exchange-correlation  potential.  The  two-component 
scheme  generates  the  ion  distribution  self-consistently 
along  with  the  electron  distribution.  The  approximations  of 
"average  atom"  models  such  as  Thomas  Fermi^  are  avoided  - 
leading  to  much  more  accurate  ionic  potential,  bound  and 
free  electron  distributions  and  the  specific  ionic 
configuration.  The  screening  effect  is  treated  in  our  model 
in  a  fully  non-linear  fashion  -  unlike  Debye-screening, 
which  is  known  to  be  inadequate  for  screening  by  bound 
electrons.  The  two-component  DFM  is  applicable  for 
arbitrary  plasma  density  and  temperature  and  is  capable  of 
treating  multielectron  ions  of  arbitrary  Z  accurately. 
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II. 


THE  METHOD  OF  CALCULATION 


a.  Isolated  Atoms  or  Multielectron  Ions 


When  the  plasma  density  is  very  low,  the  atoms  or  ions 
forming  the  plasma  can  be  considered  as  practically  non¬ 
interacting  or  isolated  from  each  other.  The  bulk  of  the 
existing  atomic  data  used  in  X-ray  laser  modeling  by 
different  research  groups  have  been  calculated  for  isolated 
ions.  In  our  modeling,  we  have  used  the  local  density 
functional  method  (DFM)  for  that  purpose.  In  actual 
calculations,  in  order  to  treat  many  electron  atoms  with 
high  z  appropriately,  relativistic  DFM  was  used.  For 
purposes  of  simplicity,  we  discuss  the  non-relativistic  DFM 
in  this  section.  The  relativistic  DFM  is  discussed  in 
Reference  4. 

The  first  part  of  the  calculation  is  to  generate  the 
energy-level  spectrum  and  the  wave-functions  of  the 
particular  atom  or  ion  of  specific  conf iguraiton .  in  the 
non-relativistic  DFM,  the  following  set  of  equations  were 
solved  iteratively  to  self-consistency: 


V  <r)  = 

a 

and 


2-  + 
r 


n(  r» 


I  dr ' 


r-r  • 


Vxc(?> 


(1) 

(2) 
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(3) 


n(r)  -  f  f • |Y.(r) |2 

In  above  n(r)  is  the  electronic  charge  density  of  the 

atom  or  multielectron  ion,  f^'s  are  the  integral  occupation 

factors  corresponding  to  the  number  of  electrons  in  each 

state  Y^(r)  with  corresponding  energy  eigenvalue  e^.  The 

atomic  potential  V_(r)  contains,  in  addition  to  the  nuclear 

and  the  electrostatic  Hartree  term,  a  contribution  arising 

from  the  exchange-corelation  effects  (V  )  due  to 

interaction  among  bound  electrons.  Let  us  note  that  the 

integer  occupation  factors  f^'s  for  the  giren  configuration 

distinguishes  this  model  from  the  "average  atom  model"  in 

which  the  occupation  factors  are  taken  to  be  those  given  by 

the  statistical  fermi  distribution  function. 

The  exchange-correlation  potential  Vxc  has  to  be 

calculated  separately  from  many-body  theory  and  used  as  an 

input  here.  In  our  calculation,  Gunnarsson-Lundquist  (G-L) 

form  for  V  was  used.  It  is  known  that  reliable  atomic 
xc 

data  is  obtained  from  G-L  exchange-correlation.  Eqs.  (1)- 
(3)  were  solved  numerically  to  self-consistency  to  obtain 
the  wavefunctions  ' s ,  the  binding  energies  e^'s  of  each 
orbital,  the  atomic  charge  density  p(r)  and  the  self- 
consistent  potential  V  (r). 

To  calculate  the  photoabsorption  by  the  atom  (or  ion), 

now  consider  the  effect  of  an  incident  time-varying 

radiation  field  E(t)-E  elwt  on  the  atom.  It  induces  a  time 

o 

dependent  atomic  density  deviation,  &n(r,t),  causing  a  time- 


7 


dependent  polarization  effect.  For  the  linear  response 
method  used  here,  it  is  convenient  to  work  with  the  Fourier 
transform: 


6n(r,t)  -  jjf  J_  Sn(r,w)e“iwt  dw  (4) 

1 

the  net  induced  density  due  to  the  external  plus  the  induced 
potential  is: 


Snind(r,w)‘J  X(r,r' ,«) [vext(r,w)+vind(r,«)]dr' ,  (5) 


where  the  induced  potential  is  given  by 


vind(t.»). 


d?.+  av»c(n(rl> 


r-r ' 


3n(  r ) 


(6) 


The  response  function  is  given  by 


X(r,r',w)-f  fiYi*(  r  )Yi(  r' )  G(r,r', 


Z  f.Y.(r)Y*(r' )G*(r,r',ei-Ro)) 


(7) 


and  thus  involves  the  wave  functions  and  energy  levels  of 
the  atom.  The  Green's  functions  are  solutions  of  the 
inhomogeneous  Schrodinger  equation  and  are  constructed 
numerically. 
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The  frequency  dependent  polarizability  o(w)  is  the 
ratio  of  the  induced  dipole  moment  to  the  external  field: 


«(«)  -  -  |  J  Z  in(r,w)dr  (8) 

Note  that  <*(«),  like  &n(r,w)  is  complex.  The  induced 
density  deviation  (  and  the  corresponding  induced  potential) 
can  have  a  phase  difference  with  respect  to  that  of  the 
applied  external  field.  Once  <x(«)  is  determined,  the 
photoabsorption  cross-section  a(«)  of  the  atom  or  the  ion  is 
calculated  from: 


a(w)  *  Im  a(w) 

c 

b.  Partial  Cross-Section 


(9) 


Eq.  (9)  gives  the  total  photoabsorption  cross-section 
which  is  the  sum  of  partial  (individual)  photoionization 
cross-sections.  To  see  the  connection  with  the  iPA-model, 
consider  the  partial  cross-section  due  to  photoionization 
from  a  specific  bound  state  ^ ( r ) . 

The  initial  atomic  state  is  represented  as 

.  Unl(r) 

Tj(  r )  -  Yx(r)  (10) 
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and  the  final  continuum  state  with  wave  vector  $  and  energy 
c  as 


▼£r) 


4n  I  A, , i 

r 


.1'  ftl'(r) 


(K)  Y1  (r) 


(11) 


The  complex  coefficients  A^'s  are  found  by  requiring  Y^.(r) 
to  behave  aasumptotically  as  an  incident  plane  wave  plus  a 
spherical  wave.  Then  the  partial  cross-section  on^  is  shown 
to  be 


2 

an^(«)  •  2  ( 21+1 )  ahco  ^le  afi 

(12) 

x  I  lAr|2  |<1  100  |l'o>l2  |  J  Pel,(r)VSCF(r,w)unl(r)dr  |2 

where  <1  100  |l'o>  is  a  clebsch-Gordon  coefficient. 

In  (12),  V  (r,w)  is  a  frequency  dependent  complex 
self-consistent  potential  generated  in  solving  iteratively 
Eq.  (5)  -  (7).  Note  that,  if  VSCF(r,w)  is  replaced  by  the 
usual  dipole  moment  operator,  one  obtains  the  conventional 
or  independent  particle  approximation  (IPA)  result.  Also, 
the  real  and  imaginary  parts  of  the  self-consistent  field 
contributes  to  the  partial  cross-section  without 
interference . 
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c.  Radiative  Decay 

The  Einstein  transition  probability  rate  is  defined  as 
the  total  probability  per  unit  time  of  an  atom  in  a  specific 
state  i  making  a  transition  to  a  final  state  j.  The  wave 
functions  and  energy  levels  for  the  atom  (or  ion)  were  first 
generated  using  the  atomic  structure  code.  The  radiative 
decay  rate  were  then  calculated  by  evaluating  the  relevant 
dipole  matrix  element  numerically  for  the  particular 
transition.  The  actual  equations  are  given  in  Ref.  12.  We 
used  (j-j)  coupling  in  our  calculations. 

d.  Radiative  Recombination 


This  is  effectively  the  inverse  of  photoionization 
process.  Thus,  once  the  photoionization  cross-sections  have 
been  computed  as  discussed  earlier,  the  readiative  capture 
cross-sections  (and  rates)  can  then  be  calculated  from  the 
process  of  detailed  balance.  The  details  can  be  found  in 
Ref.  12. 

e .  Inclusion  of  the  Effect  of  Plasma  Density  and 
Temperature . 

In  previous  sections,  we  discussed  the  calculation  of 
atomic  structure  and  processes  for  an  isolated  atom  or  ion  - 
that  is,  without  the  effect  of  the  plasma  environment.  This 
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is  valid  for  very  low  density  plasmas.  As  the  plasma 
density  increases,  several  additional  effects  due  to 
screening  shift6  of  energy  levels  (or  lowering  of  ionization 
potential)  and  modification  of  wavefunctions  and  potentials 
for  the  ion  embedded  in  high  density  plasmas  have  to  be 
considered.  We  utilized  the  two-compouent  finite 
temperature  DFM  to  incorporate  the  effects  of  the  plasma 
environment.  Consider  an  ion  (li-like  Carbon,  for  example) 
embedded  in  a  plasma  in  local  thermodynamic  equilibrium  at  a 
temperature  T.  The  charged  ion  attracts  the  free  electrons 
in  the  plasma,  while  pushing  other  charged  ions  away.  Let 
p(r)  and  n(r)  be  the  inhomogeneous  ion  and  electron 
distribution  in  the  plasma  and  p  and  n  the  corresponding 
mean  densities.  Over  a  sphere  of  radius  R  (the  correlation 
sphere)  surrounding  the  central  ion,  the  electrical 
neutrality  condition  gives 


z  f  p(r)  dr  -  z  -  [  n(r)  dr  (13) 

JR  JR 

where  z=z  -  n^,  where  n^  is  the  number  of  bound  electrons  of 
the  ion.  From  first  principles,  the  following  set  of 
coupled  equations  can  be  derived  for  electrons  and  ions  in 
the  plasma  using  the  finite  tempertaure  two-component  DFM: 


[-  \  V2  +  Ve(r)]»i(r)  -  ti*i(z) 


(14) 


12 


with 


Ve(r) 


-IW 


Z  p 

r  )-n( 

r') 

4 

r-r ' 

+  vxc(?l  -  vxc<^' 


(15) 


OD  00 

n(r)  -  Z  fb|*b(r)  |2+^2(21+1)J  dk  k2  f  ( k  )  |  *kl  (  r  )  |  2  (16) 


In  most  plasma  conditions,  the  ion  distribution  can  be  well 
represented  by 


t  ~  -V  /KT 

p( r )  *  p  e 


with  VA(r) 


[z  f  z  ptz )-n?r' )  3r'l 

Lr  J  I?-?' I  J 


(17) 

(18) 


In  Eqs.  (14)  -  (18),  we  note  that  the  effective  electron 
potential  Ve(r)  is  influenced  both  by  the  electron  and  the 
ion  charge  distributions:  n(r)  and  p(r).  The  ion  effective 
potential  Vx(r)  is  also  influenced  the  same  way.  The 
electron  density  distribution  in  eq.  (16)  is  determined  from 
the  contribution  from  the  bound  states  *b's  as  also  from  the 
continuum  states  *jci(r)  weighted  by  the  integral  occupation 
factors  fb's  for  bound  electrons  and  the  fermi  distribution 
function  f^  for  free  electrons. 
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The  model  outlined  above  incorporates  several 
improvements  over  frequently  used  "average  atom"  models  such 
as  Thermi-Fermi .  The  present  model  is  valid  at  any  plasma 
density  and  temperature,  i.e.,  both  for  weakly  and  strongly- 
coupled  plasma  conditions.  The  screening  effect  of  both 
bound  and  free  electrons  are  treated  accurately  in  this 
model,  unlike  Thomas-Fermi .  The  many-body  effects  such  as 
electron  exchange  and  correlation  are  treated  adequately. 

The  electron  degeneracy  effect  is  included  and  realistic  ion 
distribution  is  folded  into  the  calculation  and  is  made 
self-consistent  to  the  electron  distribution. 

Debye-screened  potentials  have  been  used  often  to 
describe  the  plasma  screening  effects  on  atomic  processes, 
such  as  photoionization.  For  comparison  purposes,  we 
carried  out  similar  calculations  for  a  few  representative 
density  and  temperature  using  a  simplified  model.  In  this 
model,  the  long  range  part  of  the  effective  ionic  potential 
was  represented  in  the  Debye-screened  form,  whereas  the 
inner  part  of  the  ionic  potential  was  calculated  using  the 
bound  wavefunctions.  Thus,  the  summation  over  a  large 
number  of  continuum  states  and  the  corresponding  iteration 
could  be  avoided  and  resulted  in  computational  simplicity. 
The  results  are  compared  with  the  self-consistent  DFM 
results.  Let,  us  emphasize  that  the  finite  tenperature  two- 
component  DFM  goes  far  beyond  the  Debye-screened  model. 


14 


Debye  screening  is  known  to  be  inadequate  in  a  wide  range  of 
plasma  conditions  whereas  the  DFM  treats  the  screening 
effects  in  a  fully  non-linear  self-consistent  fashion  and  is 
therefore  much  more  realistic  and  useful. 
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III.  RESULTS 


To  see  the  usefulness  of  the  time-dependent  density 
functional  method  (TD-DFM)  for  generating  photionization 
cross-sections,  we  first  discuss  the  results  for  a  few 
neutral  atoms.  In  Fig.  1,  results  of  computed 
photoionization  cross-section  for  neutral  Xenon  is  plotted 
as  a  function  of  photon  energy  near  the  4d-threshold .  The 
drastic  difference  between  the  TD-DFM( curve  A)  and  the 
conventional  IPA-model  (Curve  B)  is  clear.  The  IPA-model 

O 

does  not  reproduce  the  experimental  values  at  all  whereas 
the  present  TD-DFM  model  agrees  very  well  with  the 
experimental  data  over  this  entire  range  of  photon  energy 
from  5  to  10  Ryd.,  including  the  peak  at  about  7  Ryd.  In 
that  energy  range,  the  IPA-cross  section  shows  a  rapid 
decrease  -  in  contrast  to  the  experimental  data. 

The  primary  physical  reason  for  this  difference  in  the 
two  models  arise  from  the  polarization  effect  of  the  atom 
subject  to  the  incident  radiation.  This  collective  effect 
is  missing  in  the  iPA-calculation .  A  large  number  of 
electrons  forming  the  atom  participate  in  the  polarization 
process.  The  external  field  is  screened  in  the  energy  range 
5-6  Ryd  (and  again  in  the  range  8.5  -  9.5  Ryd.),  and  is 
antiscreened  in  the  intermediate  6-8  Ryd.  range.  The 
antiscreening  effect  produces  a  stronger  effective  field  for 
the  4d-electron  to  photoionize,  thereby  enhancing  the  cross- 
section  in  the  intermediate  energy  range  as  seen  in  Fig.  1. 
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The  3S-partial  cross-section  for  Argon  is  shown  in 

g 

Fig.  2.  The  experimental  data  in  the  range  30-65  eV  are 
depicted  by  circles  with  error  bars.  The  conventional  IPA- 
calculation  (curve  A)  again  does  not  show  the  experimentally 
seen  variation  at  all.  The  TD-DFM  model  reproduces  the 
observed  variation  including  the  experimental  Cooper 
minimum1  at  about  43  eV  and  is  in  good  agreement  with  the 
measured  cross-sections.  The  occurrence  of  the  Cooper 
minimum  is  known  to  be  due  to  the  vanishing  of  the  matrix 
element  between  bound  and  continuum  states  at  that  photon 
energy.  In  the  TD-DFM  model,  the  physical  reason  for  the 
minimum  is  that  the  induced  potential  almost  exactly  cancels 
out  the  external  field,  reducing  the  effective  field  to 
almost  zero. 

The  total  cross-section  for  Argon  atom  as  a  function 
of  the  photon  energy  is  plotted  in  Fig.  3,  and  compared  with 
experimental  data1*.  Comparison  is  also  made  of  results 
obtained  from  Hartree-Fock  (HF)  length  and  velocity 
approximations1®.  The  present  TD-DFM  model  again  best 
reproduces  the  experimental  data.  The  HF-velocity 
approximation  gives,  for  example,  a  cross-section  twice  as 
large  at  40  eV  whereas  the  HF-length  result  is  six  times 
larger  than  the  experimental  data  at  the  same  photon  energy. 

Good  agreement  with  available  experimental  data  is 
also  obtained  for  Neon  and  Krypton  by  using  the  present  TD- 
DFM  model.  This  suggests  that  the  TD-DFM  model  could  be  a 
reliable  tool  to  generate  photoionization  and 
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photoabsorption  data  for  multielectron  ions  of  different  Z. 
We,  therefore,  carried  out  the  calculations  for  several  ions 
of  experimental  interest  for  x-ray  laser  systems:  H-like  C, 
Li-like  C,  Li-like  Ne ,  Ne-like  Ar,  Ne-like  Se  etc.  For  most 
of  these  ions,  no  experimental  data  is  available  and  thus 
the  data  have  to  be  provided  by  theoretical  calculations. 

In  earlier  discussion,  it  was  pointed  out  that  the 
collective  effects  of  a  large  number  of  bound  electrons  is 
reponsible  for  the  substantial  difference  between  TD-DFM  and 
Independent  Particle  Approximation  (IPA)  results.  For 
highly  stripped  ions  (such  as  the  ones  listed  above), 
however,  both  models  tend  to  produce  similar  results.  For 
H-like  C,  for  example,  there  is  only  1  bound  electron.  The 
TD-DFM  and  the  IPA-calculations  therefore  give  the  same 
result  (Fig.  4).  For  Li-like  Carbon  with  3  bound  electrons, 
the  photoionization  cross-sections  are  slightly  different 
for  the  photon  energy  range  370-600  eV  (Fig.  5).  Li-like 
Neon  also  has  three  bound  electrons.  However,  because  of 
higher  nuclear  charge  (Z«10),  those  electrons  are  tightly 
bound  (more  tightly  than  in  neural  Neon)  and  thus  do  not 
participate  as  effectively  in  the  polarization  process. 

Thus,  the  TD-DFM  and  IPA  calculations  lead  to  essentially 
the  same  results  for  Li-like  Ne  (Fig.  6(a)  and  6(b))  and 
Ne-like  Se  (Fig.  8).  For  Ne-like  Argon,  the  polarization 
process  is  somewhat  more  effective  and  leads  to  results  that 
are  somewhat  different  from  IPA  results  over  the  photon 
energy  range  450-600  eV  (Fig  7).  To  see  the  variation  of 


18 


the  contribution  of  each  electronic  state  to  the  total 
cross-section,  the  computed  values  for  Li-like  C,  Li-like  Ne 
and  Ne-like  Ar  are  given  in  Tables  i-m. 

The  radiative  decay  rate  of  ions  of  specific 

configuration  is  another  kind  of  atomic  data  necessary  for 

atomic  modeling  for  X-ray  lasers.  Once  the  wavef unctions 

and  energy  levels  of  the  given  multielectron  ion  has  been 

generated  from  the  atomic  structure  code,  the  radiative 

decay  rate  (or  Einstein  transition  probability)  can  be 

12 

calculated  by  evaluating  the  relevant  matrix  element  .  In 
Table  IV,  the  computed  decay  rates  of  some  of  the 
transitions  for  Neon-like  Krypton  is  presented.  (j-j) 
coupling  has  been  used  in  these  calculations.  The  ground 
state  of  Ne-like  Krypton  is  lS22S22p6  with  total  J=0.  All 
energies  in  Table  IV  are  given  relate-  to  the  ground  state 
energy  of  the  ion  which  is  taken  to  be  Zero.  We  note  from 
the  table  that  the  radiative  decay  rates  vary  substantially 
depending  on  the  particular  excited  ionic  state  and  the 
transition  considered. 

Next,  we  incorporated  the  effect  of  plasma  density  and 
temperature  into  our  model  and  investigated  the  resulting 
modification  of  photoionization  cross-sections  of  ions  in 
plasmas.  The  plasma  effects  were  studied  first  by  using  an 
approximate  model  in  which  the  long  range  part  of  the  ionic 
potential  was  taken  in  the  form  of  Debye-screened  potential 
(for  the  given  plasma  density  and  temperature).  The  inner 
part  of  the  potential  was  constructed  numerically  using  the 
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bound  state  wavef unctions  as  in  Eq.  (15).  This  avoids  the 
large  number  of  iterations  with  the  continuum  wavefunctions 
and  thus  leads  to  computational  simplicity* ^ . 

We  applied  this  simplified  model  to  H-like  and  Li-like 

Carbon  ions  in  plasmas.  The  results  are  shown  in  Figures 

9-12.  The  plasma  densities  considered  are  in  the  range 
22  23  3 

10  -  10  electrons/cm  and  temperatures  are  in  the  range 
30  -  50  eV.  The  corresponding  Debye  screening  lengths  are 
in  the  range  of  2  -  5  au. 

From  Fig.  9,  we  see  by  comparing  with  isolated  ion 
cross-sections,  that  the  threshold  'for  photoionization  from 
ls-level  of  H-like  C  shifts  substantially  (from  18  eV  to 
about  15.2  eV)  for  Debye  length  XD*=2  au.  This  is  important 
in  the  context  of  calculation  of  photoionization  rates, 
which  is  calculated  by  integrating  the  cross-section  weighed 
by  the  electron  distribution  function  (usually  Maxwell- 
Boltzmann  distribution)  over  the  entire  energy  range.  The 
calculated  rates  in  the  two  cases  would  therefore  be  very 
different  -  as  seen  from  Fig.  9. 

In  Fig.  10,  the  corresponding  results  of  2S 
photoionization  cross-section  for  H-like  C  is  plotted.  Fig. 

11  shows  the  3d  partial  cross-sections  of  the  same  ion  at 
two  different  plasma  conditions  (XD=5  and  2  au)  and  are 
compared  with  the  isolated  atom  cross-section  (curve  A). 

The  increasing  shift  of  the  ionization  threshold  with 
decreasing  XD  (i.e.,  with  increasing  density)  is  clear  from 
the  Figure  11.  For  Li-like  Carbon,  the  2S  partial  cross- 
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section  at  two  different  plasma  conditions  were  computed 
using  the  simplified  model  and  are  shown  in  Fig.  12. 
Comparison  with  isolated  ion  cross-section  shows  substantial 
modification  of  photoionization  cross-section  at  these 
densities  and  temperatures  and  significant  shift  of 
ionization  threshold  due  to  the  effect  of  the  surrounding 
plasma . 

The  full  self-consistent  calculation  using  the  two- 
component  DFM-model  was  also  performed  for  Li-like  Carbon  in 
two  different  plasma  conditions.  The  results  are  shown  in 
Fig.  13  (  Xd«5  au).  and  Fig.  14  ( XD«=2  au).  Curve  B  in  both 
figures  correspond  to  data  obtained  from  the  simplified 
model  whereas  curve  C  is  the  result  of  full  self-consistent 
DFM  calculation.  The  DFM  model  generates  more  accurate 
shifted  energy  levels;  thus  the  shift  of  the  threshold  of  2S 
cross-section  relative  to  the  isolated  ion  threshold  is  more 
realistic  in  this  model.  In  addition,  the  self-consistent 
DFM-model  is  seen  to  generate  photoionization  cross-sections 
that  are  substantially  different  from  simplified  model.  The 
calculated  rates  resulting  from  these  cross-sections  also 
differ  substantially. 
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IV. 


SUMMARY 


Some  of  the  atomic  processes  important  for  overall 
modeling  of  x-ray  lasers  have  been  investigated  in  this 
work,  performed  over  the  period  June  1986  through  July  1987. 
The  atomic  structure  code  for  multielectron  ions  of 
arbitrary  Z  is  of  central  importance  in  these  computations. 
The  code  we  have  used  for  this  purpose  is  based  on  local 
density  functional  method  and  provides  a  fast,  efficient, 
accurate  tool  to  generate  energy  levels,  wave-functions, 
etc.  for  multielectron  ions  and  atoms.  Proper  relativistic 
treatment  is  included  in  this  code  for  high  Z  ions.  We 
carried  out  computations  of  photoionization  and 
photoexcitation  cross-sections,  radiative  recombination, 
radiative  decay  rates  for  a  variety  of  ions  of  different  Z 
over  a  wide  range  of  parameters  of  interest  to  X-ray  lasers. 

It  is  demonstrated  that  the  time-dependent  linear 
response  method  within  the  framework  of  density  functional 
theory  can  provide  reliable  photoprocess  (and  recombination) 
data  for  various  atoms  and  ions  of  experimental  interest. 
This  model  is  particularly  useful  in  those  situations  where 
conventional  independent  particle  models  fail  to  provide 
accurate  data.  The  mechanism  of  time-dependent  polarization 
of  the  atom  is  seen  to  be  important  in  describing  the 
experimental  results.  As  a  practical  point,  the  computer 
code  based  on  the  time-dependent  model  is  fast  and 
efficient,  capable  of  generating  a  large  number  of  data  in  a 
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short  time  (for  example,  cross-sections  for  10  photon 
energies  for  a  medium-z  ion  takes  about  3  mins,  of  C.p.u. 
time  on  a  Cray-XMP  computer). 

We  also  studied  the  effects  of  plasma  density  and 
temperature  on  atomic  structure  and  processes  such  as 
photoionization  and  spontaneous  transitions.  The  model  and 
the  computer  code  developed  for  this  purpose  is  capable  of 
treating  screening  effects  in  a  fully  non-linear  fashion  and 
thus  goes  beyond  the  often  used  Debye-screening  model.  The 
model  represents  a  significant  improvement  over  the  "Average 
Atom"  models  (AAM)  in  a  number  of  ways.  The  specific  ionic 
configuration  (with  integral  occupation)  is  maintained  in 
our  model-unlike  AAM's.  Both  the  bound  and  free  electrons 
are  treated  fully  quantum  mechanically  in  this  model  (in 
contrast  to  AAM's  such  as  Thomas-Fermi )  and  therefore,  lead 
to  more  accurate  electronic  charge  distribution  and 
effective  ionic  potentials.  The  ion  distribution  in  our 
model  is  treated  self-consistently  with  the  electron 
distribution  via  coupled  equations.  The  present  model  is 
applicable  in  arbitrary  plasma  conditions  and  is  therefore, 
useful  in  a  wide  range  of  densities  and  temperatures  where 
models  such  as  Thomas-Fermi  or  Debye-Huckel  are  known  to  be 
inadequate.  The  calculations  done  here  show  that,  with 
increasing  plasma  density,  the  atomic  data  (such  as 
photoionization  cross-section)  can  be  substantially  modified 
from  the  isolated  atomic  values.  The  effects  of  plasma 


density  and  temperature  thus  need  to  be  included  for 
realistic  atomic  modeling.  The  models  presented  here 
achieve  that  objective  systematically. 

The  work  performed  have  thus  provided  improved  atomic 
modeling  capability  for  X-ray  laser  plasmas  and  should  be 
useful  for  application  to  different  lasing  schemes  such  as 
photopumping  or  radiative  recombination.  In  addition,  the 
models  and  the  computer  codes  are  useful  tools  for  studies 
of  radiative  properties  of  hot  plasmas  of  different 
densities,  opacity  calculations  and  plasma  diagnostics. 
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TABLE  I 


Li-Like  C 


Photon 
Energy  (ev) 


Electronic 

States 


Photo ionization 
_ ( Megab 


250 

2S  1/2 

I  PA 
0.0424 

300 

2S  1/2 

0.0280 

IS  1/2 

0.5968 

350 

2S  1/2 

0.0193 

Total 

0.6161 

IS  1/2 

0.4354 

400 

2S  1/2 

0.0141 

Total 

0.4495 

IS  1/2 

0.3088 

460 

2S  1/2 

0.0100 

Total 

0.3188 

IS  1/2 

0.2499 

500 

2S  1/2 

0.0081 

Total 

0.2S80 

IS  1/2 

0.1954 

550 

2S  1/2 

0.0065 

Total 

0.2019 

IS  1/2 

0.1556 

600 

2S  1/2 

0.0051 

Total 

0.1608 

Cross-Section 
rn) _ 


TD-DFM 

0.0370 

0.0275 

0.5969 

0.0190 

0.6159 

0.4444 

0.0141 

0.4585 

0.3203 

0.0101 

0.3304 

0.2613 

0.0083 

0.2695 

0.2059 
0.0066 
0 . 212S 

0.1650 

0.0053 

0.1703 
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TABLE  II 


Li-Like  Ne 


Photon 
Energy  (ev) 

Electronic 

States 

Photoionization  Cross-Secti 

(Meqabarn) 

I  PA 

TD-DFM 

230 

2S  1/2 

0.2108 

0.2080 

240 

2S  1/2 

0.1960 

0.1932 

260 

2S  1/2 

0.1673 

0.1649 

280 

2S  1/2 

0.1451 

0.1430 

300 

2S  1/2 

0.1265 

0.1247 

320 

2S  1/2 

0.1116 

0.1099 

340 

2S  1/2 

0.0986 

0.0971 

360 

2S  1/2 

0.0877 

0.0863 

1040 

2S  1/2 

0.0080 

0.0099 

1060 

2S  1/2 

0.0076 

0.0074 

1080 

2S  1/2 

0.0073 

0.0082 

IS  1/2 

0.0379 

0.0395 

2000 

2S  1/2 

0.0015 

0.0016 

Total 

0.0395 

0.0411 

IS  1/2 

0.0359 

0.0375 

2040 

2S  1/2 

0.0014 

0.0015 

Total 

0.0373 

0.0390 
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TABLE  III 


Ne-Like  Ar 


Photon 
ergy  (ev) 

Electronic 

States 

Photoioni zation 

(Meqab 

395 

2p  3/2 

Ip a  — 

1.0061 

400 

2p  1/2 

0.4984 

2p  3/2 

0.9905 

Total 

1.4889 

420 

2p  1/2 

0.4339 

2p  3/2 

0.8593 

Total 

1.2931 

440 

2p  1/2 

0.3952 

2p  3/2 

0.7807 

Total 

1.1759 

460 

2p  1/2 

0.3509 

2p  3/2 

0.6938 

Total 

1.0447 

480 

2S  1/2 

0.1981 

2P  1/2 

0.3188 

2p  3/2 

0.6297 

Total 

1.1466 

c00 

2S  1/2 

0.1901 

2p  1/2 

0.2867 

2p  3/2 

0.5656 

Total 

1.0424 

520 

2S  1/2 

0.1769 

2p  1/2 

0.2604 

2p  3/2 

0.5139 

Total 

0.9512 

540 

2S  1/2 

0.1682 

2p  1/2 

0.2368 

2p  3/2 

0.4668 

Total 

0.8719 

560 

2S  1/2 

0.1497 

2p  1/2 

0.1971 

2p  3/2 

0.3884 

Total 

0.7351 

600 

2S  1/2 

0.1422 

2p  1/2 

0.1804 

2p  3/2 

0.3553 

Total 

0.6779 

Cross-Section 

rn) _ 

TD-DFM 

0.9481 

0.5060 

1.0065 

1.512S 

0.4465 

0.8855 

1.3321 

0.4183 

0.8273 

1.2457 

0.3774 

0.7468 

1.1243 

0.2073 

0.3471 

0.6863 

1.2407 

0.1995 

0.3159 

0.6239 

1.1393 

0.1863 

0.2893 

0.5717 

1.0473 

0.1776 

0.2657 

0.5243 

0.9676 

0.1589 

0.2240 

0.4419 

0.8247 

0.1513 

0.2063 

0.4068 

0.7643 
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TABLE  IV 


Neon-Like  Krypton 


Initial  State 


Energy* 

(ev) 


Final  State 


Energy* 
lev)  _ 


Radiative 
Decay  Rate** 


2p53s (3/2, 1/2) J=1 

1652.94 

2p6(0,0)  J=0 

0.00 

5.04E+12 

2p53d (3/2,3/2)J=l 

1782.71 

2p6  (0, 0) J=0 

0.00 

2.58E+11 

2p53p (3/2, 1/2) J=1 

1701.40 

2p^3s 

(3/2,  1/2) J=2 

1649.75 

7.87E+09 

2p~*4p 

2282.36 

2p53s 

(3/2, 1/2) J=2 

1649.95 

6. 97E+11 

2p54  f 

2324.58 

2p^3s 

(3/2,  1/2) J=2 

1649.95 

9.24E+07 

2p54s 

2237.98 

2P53p 

(3/2,  1/2) J=1 

1701.40 

2.54E+11 

2p54d 

2310.79 

2p53p 

(3/2, 1/2) J=1 

1701.40 

3.27E+11 

2p53s (1/2, 1/2) J=1 

1707.72 

2p53p 

(3/2, 1/2) J=2 

1704.20 

2.03E+02 

2p53d(3/2, 3/2) J=1 

1781.71 

2p53p 

(3/2,  1/2) J=2 

1704.20 

2.54E+09 

2p54p 

2282.36 

2p53d 

(3/2,  3/2) J=0 

1780.34 

6.84E+09 

2p54f 

2324.58 

2p53d 

(3/2, 3/2) J=0 

1780.34 

1.36E+11 

2p54d 

2310.79 

2p^4p 

2282.36 

5.68E+09 

*A11  energies  are  measured  relative  to  the  ground  state: 
2p^(0,0)J=0.  The  numbers  in  parenthesis  are  the  j-value  of  the 
parent  ion  and  that  of  the  outer  electron. 
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Figure  1 
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Figure  3 


C9W)  N0I133S-SS0yD  lUiOi 
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SOLID  CURVE:  TIME-DEPENDENT  DENSITY  FUNCTIONAL  CALCULATION. 
X:  EXPERIMANTAL  DATA  FROM  J. A. R. SAMSON,  ADVAN.  ATOM.  MOL. 
PHYS,  2,  178(1966).  +:  HARTREE-FOCK  (VELOCITY  APPROX.)  AND 
0:  HARTREE-FOCK  (LENGTH  APPROX.)  BY  KENNEDY  &  MANSON, 

PHYS.  REV.  A  5,  227(1972) . 
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Figure  5 
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SOLID  CURVE:  INDEPENDENT  PARTICLE  MODEL. 

DASHED  CURVE:  TIME -DEPENDENT  DENSITY  FUNCTIONAL  CALCULATION. 


2S  PfiRTIRL  CROSS-SECTION 


Figure  6 (b) 
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CURVE  A:  INDEPENDENT  PARTICLE  MODEL.  CURVE  B:  TIME- 
OEPENDENT  DENSITY  FUNCTIONRL  CALCULATION. 
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Figure  7 


SOLID  CURVE:  INDEPENDENT  PARTICLE  MODBL. 

DASHED  CURVE:  TIME-DEPENDENT  DENSITY  FUNCTIONAL  CALCULATION. 
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Figure  8 
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SOLID  CURVE:  INDEPENDENT  PARTICLE  MODEL. 


DEBYE  LENGTH  =  2  ( AU ) 


2S  LEVEL  CROSS  SECTIONS 


Figure  10 
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ZERO  DENSITY,  B:  DEBYE  LENGTH  =  2  (AU) 
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Figure  11 
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ZERO  DENSITY,  B:  DEBYE  LENGTH  =  5  (AU) 

DEBYE  LENGTH  =  2  (AU). 


Figure  12 
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ZERO  DENSITY,  B:  DEBYE  LENGTH  =5  (AU) 

n  rtifn  t  min  m  *  *  n  /  ■  t  r  \  ' 


CARBON  IV 
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ZERO  DENSITY,  B:  DEBYE  LENGTH 
SELF  CONSISTENT  FOR  DEBYE  LEN< 


CARBON  IV 


Figure  14 
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ZERO  DENSITY,  B:  DEBYE  LENGTH  =  i 
SELF  CONSISTENT  FOR  DEBYE  LENGTH 
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FOR 
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INTRODUCTION 


I . 

Electron  collisional  excitation  is  one  of  the  important 
mechanisms  through  which  population  inversion  of  ionic  energy 
levels  leading  to  lasing  in  the  soft  x-ray  region  can  be  achieved. 
Several  research  groups  had  suggested  schemes  based  on  electron 
collisional  excitation  processes  as  the  pumping  mechanism1-^. 
Recent  experimental  observation  at  Livermore  of  3p-3s 
amplification  in  Ne-like  Selenium  at  206  and  209  A  using  a  laser 
driven  exploding  foil  target  demonstrates  the  usefulness  of 
electron  collisional  excitation  (as  also  the  radiative 
recombination  processes)  as  a  pumping  mechanism  leading  to  lasing. 

Aside  from  the  application  to  x-ray  lasers,  electron 
collisional  ionization  and  excitation  processes  are  among  the 
various  important  atomic  processes  in  overall  modeling  of 
radiative  properties  of  plasmas  in  a  wide  range  of  experimental 
conditions.  In  the  non  local  thermodynamic  equilibrium  (NLTE) 
conditions,  which  frequently  occur  in  plasmas,  the  detailed  rate 
equation  approach  to  modeling  requires  a  large  variety  of  atomic 
data.  Electron  collisional  ionization  and  excitation  cross- 
sections  and  rates  (along  with  other  bound-bound,  bound-free  and 
free-free  processes  such  as  photoionization,  photoexcitation, 
radiative  recombination  and  decay,  etc.)  are  necessary  inputs  to 
these  calculations.  Accurate  calculations  of  electron  collisional 
ionization  and  excitation  processes  are  therefore  needed  for  the 
different  modeling  mentioned  above.  Also,  even  though  some 
experimental  data  for  scattering  cross-sections  exist  for  certain 
atoms  and  ions;  there  is  a  lack  of  experimental  data  for  many 
different  types  of  atoms  and  ions  over  a  wide  range  of  energy. 
Thus,  the  required  data  have  to  be  provided  by  theoretical 
calcu) scions . 

In  theoretical  calculations  of  electron  scattering  cross- 
sections,  semi-classical  impact  (SCI)  approximation  has  been  used 
in  some  cases  to  generate  the  data^.  The  SCI-method,  however, 
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turns  out  to  be  inadequate  in  many  cases.  Firstly,  the  SCI-method 
requires  the  choice  of  an  impact  parameter  cutoff  for  which  there 
is  no  absolute  prescription.  This  arbitrariness  can  lead  to 
certain  errors.  The  method  itself  is  restricted  to  dipole  allowed 
excitations.  Also,  for  large  line  strengths,  the  assumption  of 
weak  coupling  in  SCI-method  is  inadequate.  Improving  the  first 
Born  approximation  requires  a  partial  wave  analysis. 

In  view  of  these,  it  is  necessary  to  employ  methods  such  as 
the  Distorted  Wave  with  exchange  (DWEX)1®  or  the  Close  Coupling 
(CC)  scheme11  to  carry  out  the  calculations  of  electron 
collisional  cross-sections.  The  relative  merits  of  these  methods 
will  be  clear  from  the  outline  of  these  schemes  presented  in  the 
next  section.  In  general,  the  close  coupling  method  can  provide 
data  of  high  accuracy;  although  in  many  cases,  the  DWEX  can 
generate  data  within  a  few  percent  of  the  close  coupling  results. 

For  this  6-month  project  (August  1987  -  February  1988) ,  we 
carried  out  computation  of  electron  impact  ionization  and 
excitation  cross-sections  and  rate  coefficients  for  different  ions 
of  interest  over  a  wide  range  of  incident  electron  energy.  Most 
of  the  calculations  were  performed  using  the  Distorted  Wave  with 
exchange  (DWEX)  method  -  since  this  provides  a  speedy  and  fairly 
accurate  tool  for  these  calculations.  The  outline  for 
computational  procedure  and  relative  merits  of  the  DWEX  and  close¬ 
coupling  methods  are  given  in  Section  II.  Some  selected  results 
for  collisional  ionization  cross-sections,  rate  coefficients  for 
different  electronic  transitions  for  ions  of  interest  are 
presented  in  Section  III.  The  present  work  forms  a  useful 
component  for  the  overall  atomic  modeling  capability  for  x-ray 
laser  and  laboratory  plasmas. 


II.  THE  METHOD  OF  CALCULATION 


In  this  section,  we  will  briefly  outline  the  Distorted  Wave 
and  the  Close-coupling  methods  of  electron-atom  (or  ion) 
scattering.  The  Close-coupling  (CC)  method  is  computationally 
more  complex  and  requires  substantial  computer  time.  Use  of  good 
target  wave  functions  are  essential  for  accuracy  in  calculation. 

A  short  discussion  of  target  wavefunctions  is  therefore,  included. 

A.  The  Close  Coupling  Method 

Consider  an  atomic  system  containing  (N+l)  electrons. 

The  electron  coordinates  are  xn  =  (rn,  an)  ,  n=l  to  (N+l)  ,  where  Fn 
is  a  position  vector  and  pn  is  a  spin  coordinate.  States  of  the  N 
electron  system  are  referred  to  as  "target"  states. 


Let  x 


D  Si  Li  "i 


be  a  target  state 


depending  on  all  coordinates  except  xn .  The  angular  momenta  of 
the  target  are  specified  by  Ms^  ML^  and  is  the  parity  of 

the  target.  Let  the  set  of  quantum  numbers  IN  Sj_  1^  be 
collectively  denoted  by  i.  We  will  refer  to  i  as  the  channel 
index.  In  the  close-coupling  approximation,  the  complete 
wavefunction  is  taken  to  be 


1 

*  =  A  Z  X, 

i  =  l 


F.  (r  ) 
1  n 


+ 


J 

Z 

j  =  l 


$  . 

3 


(1) 


where  A  is  the  antisymmetrization  operator,  $j's  are  the  square 
integrable  functions  for  the  entire  system,  the  radial  functions 
F-^(r)  and  the  coefficients  Cj's  are  fully  optimized.  The  F^(r)'s 
are  made  orthogonal  to  all  radial  functions  with  1=1^  used  to 
construct  the  states  x^  and  The  set  of  functions  <£j  must  be 

such  that  the  imposition  of  these  orthogonality  conditions  does 
not  imply  a  constraint  on  the  total  wavefunction  lP. 
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Let  E  be  the  total  energy,  Ej_  the  energy  for  the 
target  state  i  and  the  energy  of  the  added  electron.  Then 

E  =  Ei  +  (2) 

At  energies  giving  bound  states  for  the  whole  system  we  have 
<  0  for  all  i  and  the  radial  functions  F^(r)  go  to  zero 
exponentially  in  the  limit  r-*co.  The  wavefunctions  for  bound 
states  are  normalized  to  unity.  The  collision  states  are  such 
that  >  0  for  some  values  of  i.  Channel  i  is  said  to  be  open 
if  ei  >  0  and  closed  if  e ^  <0.  Let  the  E^  be  ordered  such  that 
E1  <  e2  <  *■'  <  Ei*  For  a  given  E,  we  have  IQ  open  channels  and 


*i  >  0  f0r  1  *  1  CO  :o  (3) 

<  0  for  i  =  (IQ+1)  to  I. 

For  numerical  work,  we  first  calculate  real  functions 
F.^, (r)  satisfying  the  boundary  conditions 

Fii' (r)  ~  si (r) 

r-»oo 

F. ., (r)  ~  0  for  i  =  (I  +1)  to  I, 

r-*0 

with  i'=l  to  IQ.  The  functions  S^(r)  and  (r)  are  Coulomb 
functions  and  K  is  the  reactance  matrix. 

Further,  construct  the  Coulomb  functions 

=  c  ±  is  (5) 


5  (i,  i'  )  +  c i  (r)  K  ( i ,  i  '  )  , 
for  i  =  1  to  I 


(4) 


with  the  asymptotic  forms 
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exp  (±i$) 


(6) 


with 


<f>“~  (jrk) 


-1/2 


$  =  kr  -  \  Iff  +  (!)  In  (2kr)  +  arg  r(l+l  -  iz/k)  (7) 

where  z  is  the  charge  on  the  target  system,  z  =  (Z-N)  with  Z  as 
the  nuclear  charge. 

The  cross-sections  can  be  determined  from  S-matrix 
which  is  related  to  asymptotically  as: 


The  collision  strengths  are  usually  expressed  in  terms  of  the 
square  of  the  appropriate  T-matrix  elements  and  therefore  to  the 
S-matrix  elements. 


B .  Distorted  Wave  Method 

It  is  well  known  that  the  DW  approximation  is 
appropriate  whenever  there  are  a  sufficiently  large  number  of 
angular  momentum  states  contributing  to  the  cross-section.  For 
highly  stripped  ions,  this  method  can  be  used  for  all  angular 
momentum  contributions. 

Let  be  the  trial  wavefunction  of  the  total  system  with 
the  asymptotic  form 


5 


-»  kp1/2$r  r j  sin(Xj+  Tp) 

-1/9  - 1 

+  I  p(f,r' )  kr,  $p,  r j  Cos(x^+  Tr,  ) 


Lj  T' 


Where  kpZ  is  the  energy  in  Rydbergs  of  the  incident 
electron  in  the  channel  r=«LaSalQ,LTST  with  ±a  the  orbital  angular 
momentum  of  the  colliding  electron,  LT  and  ST  are  orbital  and 
spin  angular  momenta  of  the  total  system  (atom  +  colliding 
electron) ,  and  Sa  the  orbital  and  spin  angular  momentum  of  the 

colliding  electron.  p(r,r')  are  the  elements  of  the  p-matrix 

which  in  turn  are  related  to  S-matrix  of  scattering  theory.  The 
matrix  element  p(r,T' )  of  the  trial  matrix  pfc  is  given  by 


pd\r')  =  PL(r,r' 


-  2  <  «fp  I  H-E  I  4>p,  > 


where  H  is  the  Hamiltonian  of  the  total  system  and  E  is  the  totaj 
energy . 

The  trial  wavefunctions  ’I'p  are  constructed  from  one- 
electron  orbital  wavefunctions  of  the  form 


*nlm  =  hm10'^  Pnllrl/r  (11> 

for  atomic  electrons  and 

*klm  =  hm'®'*1  Fkl<rl/r  <12> 

for  the  colliding  electron.  Radial  functions  Pn^  should  satisfy 
the  corresponding  Hartree-Fock  equations  for  atomic  states  aLaSa 
or  PLpSp.  F^  is  a  solution  of  the  equation 


d2  1(1+1)  2Z  ] 

- *■  -  - 5 -  + - 2  V  (r )  +  kr^  F,  .  =  0 

dr  r  r  l  kl 
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with  the  asymptotic  form 


Fkl  kr~1/2  sin(*  +  rp)  (14) 

where  Zn  is  the  nuclear  charge  and  V(r)  is  given  by 

V(r)  =  i  f  I  P  ■ 2  (y)  dy  +  f  J-  I  P.2(y)dy  (15) 

J  0  i  JrYi  1 

with  p.» s  the  radial  functions  of  all  atomic  electrons  associated 
with  the  channel  T.  The  orthogonality  of  radial  functions  are 
incorporated  in  the  procedure  via  the  method  of  Lagrange 
multipliers.  The  atomic  potential  in  eq. (15)  is  supplemented  by 
the  exchange  and  correlation  potential  as  in  the  Slater's 

1  o 

xar-method  .  In  our  calculation,  the  one-electron  orbital 
wavefunctions  are  generated  self-consistently  by  iterative 
procedure  by  using  the  atomic  structure  code  developed  in  our 
earlier  work12. 

C .  The  Target  Functions 


Use  of  good  target  wavefunctions  are  essential  for 
accurate  scattering  calculations.  Radial  functions  of  two  types 
can  be  employed;  spectroscopic  orbitals,  nl,  of  the  type  which 
one  would  have  in  a  simple  independent  particle  model  and 

correlation  orbitals  nl,  which  are  included  to  give  better 
accuracy . 


In  close  coupling  calculations  it  is  convenient  to  use 
the  same  orthonormal  sets  of  orbitals  for  all  of  the  functions  Xj, 
and  $j.  For  illustration,  consider  the  simple  case  of  He-like 
targets  for  which  one  may  wish  to  include  target  states: 


Is  2p  1P°  and  Is  2p  3P°. 
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Complete  optimization  of  the  orbitals  using  the  Hartree-Fock 
method  will  give  two  radial  functions:  (2p)1,  for  1P(-)  and  (2p)  3 

O  A  ... 

for  P  ,  which  may  be  significantly  different.  One  procedure  is 
to  use  one  spectroscopic  function  2p  and  one  correlation  function 

3p  and  to  take  the  target  functions  to  be  linear  combinations  of 

the  two  configurations  Is  2p  and  Is  3p.  With  a  suitable  choice 

O  — 

of  the  P  and  3p  functions,  this  can  be  equivalent  to  the  use  of 
the  fully  optimized  functions  (2p)  ^  and  (2p>3.  Further 
configurations  can  be  added  if  better  accuracy  is  required.  The 
radial  functions  are  calculated  numerically  from  the  Schrodinger 
equation  with  an  effective  atomic  potential.  The  effective 
atomic  potential  in  our  calculations  consists  of  the  nuclear 
term,  the  electrostatic  Hartree  term  plus  the  exchange- 
correlation  potential  that  accounts  for  the  interaction  of  the 
many  electrons  forming  the  atom  or  ion. 
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III.  RESULTS 


The  Distorted  Wave  with  Exchange  method  was  employed  to 
compute  the  electron  collision  strengths  and  cross-sections  for 
Li-like  ions  of  interest  to  x-ray  laser  plasmas.  The 
calculations  were  carried  out  over  a  large  range  of  incident 
electron  energy  starting  from  near  the  threshold  of  the 
particular  electronic  transition  considered. 

Individual  contributions  of  the  partial  waves  are  listed  in 
Table  1  to  illustrate  the  variation  of  partial  scattering  cross- 
section  with  increasing  angular  momentum.  Examination  of  Table  1 
shows  that  for  each  sec  of  energy,  the  partial  cross-section 
first  increases  with  increasing  1,  peaks  at  an  intermediate  1 
value  and  falls  off  gradually  at  higher  1.  The  total  cross- 
section  has  the  highest  value  near  threshold  and  decreases  slowly 
with  increasing  electron  energy.  This  feature  is  also  seen  from 
variation  of  collisional  excitation  cross-sections  with  incident 
electron  energy  for  the  other  transitions  of  Li-like  Al  as  shown 
in  Tables  2  and  3. 

Electron  collisional  excitation  rate  coefficients  were  also 
calculated  for  different  transitions  over  a  range  of  temperature 
appropriate  for  x-ray  laser  plasmas  and  typical  z-pinch  plasmas. 
The  excitation  rates  were  obtained  by  integrating  the  calculated 
cross-sections  weighted  by  the  Maxwell-Boltzmann  distribution 
over  the  entire  energy  range.  As  seen  from  Tables  4,  5  and  6, 
the  excitation  rate  coefficients  initially  increase  rapidly  with 
temperature,  goes  through  a  peak  and  decreases  slowly  with 
increasing  temperature.  The  maximum  of  excitation  rates 
indicates  the  temperature  at  which  electron  collisional 
excitation  can  substantially  populate  the  upper  levels  leading  to 
population  inversion  and  lasing. 

Calculations  for  other  ions  of  lithium  isoelectronic 
sequence  such  as  CIV  have  also  been  done.  We  will  briefly 
discuss  the  main  aspects  of  those  results.  Burke  et  al^  were 
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the  first  to  apply  the  close-coupling  method  to  lithium  like 
ions.  For  CIV,  calculation  of  2s-»2p  collision  strengths  show 
satisfactory  agreement  between  theory  and  experiment. 

Exchange  effects  turn  out  to  be  important  for  CIV  from 
comparison  with  experimental  data.  Distorted  Wave  without 
exchange  calculations  give  results  that  are  at  least  15%  higher 
than  experimental  data.  2-state  close-coupling  calculations,  on 
the  other  hand,  agree  with  the  experimental  data  very  well. 

The  effect  of  coupling  to  higher  states  for  2s-»2p  collision 
strength  calculation  is  much  less  in  CIV  compared  to,  for  example 
Li-like  Be.  Since  the  energy  difference  between  2s  and  2p  states 
increases  linearly  with  Z  and  the  energy  difference  between  2s 

,  o 

and  higher  n  states  increases  approximately  as  Z  ,  the  coupling 
becomes  less  important  as  Z  increases.  The  average  effect  of 
resonance  contribution  to  collision  strength  is  an  enhancement  of 
3%  in  the  energy  region  below  the  n=3  levels.  Resonances  are 
found  to  have  the  most  significant  effect  on  D-wave  (1=2) 
contribution  which  are  small  compared  with  the  total  collision 
strength . 
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IV.  CONCLUSIONS 


Accurate  calculations  of  electron  collisional  excitation 
and  ionization  cross-sections  and  rate  coefficients  of  different 
atoms  and  ions  are  important  input  to  overall  atomic  modeling  of 
radiative  properties  of  plasmas.  It  is  particularly  useful  in 
the  context  of  collisional  excitation  scheme  for  x-ray  lasers. 

In  most  investigations  to  achieve  lasing  with  ionized  plasmas  of 
different  z,  available  experimental  data  for  collisional  cross- 
sections  and  rates  are  limited.  Thus,  the  necessary  data  have  to 
be  provided  by  theoretical  calculations.  The  Distorted  Wave  with 
exchange  method  is  seen  to  be  a  reliable  and  speedy  technique  for 
the  purpose.  In  specific  cases,  such  as  near  the  threshold  for 
low  Z  elements,  close-coupling  calculations  can  provide  data  of 
high  accuracy.  The  present  work  on  collisional  processes,  (along 
with  our  previous  work  on  photoionization,  photoexcitat ion, 
radiative  decay  and  recombination  processes)  forms  a  useful 
component  for  the  overall  atomic  modeling  capability  for  x-ray 
laser  and  laboratory  plasmas. 

At  very  low  plasma  density,  collisional  excitation  and 

deexcitation  rates  are  small,  radiative  decay  is  then  a  dominant 

process.  As  the  plasma  density  increases,  the  collisional  rates 

increase  rapidly  and  collisional  radiative  equilibrium  can  exist. 

These  in  turn,  have  their  effects  on  the  relative  abundance  of 

ionic  populations  as  well  as  on  the  population  inversion  of 

levels  and  thus  affect  lasing.  The  modification  of  collisional 

excitation  and  ionization  cross-sections  etc.  with  increasing 

plasma  density,  therefore,  is  an  additional  problem  area  that 

requires  further  investigation.  In  the  present  6-month  project; 

this  aspect  has  not  been  addressed.  Preliminary  calculations, 

using  screened  atomic  potential  to  include  the  plasma  effects, 

indicate,  however,  that  for  typical  x-ray  laser  plasma  condition 
.  o  n  o  1  _  o 

(with  electron  density  10  -  10  cm  and  for  highly  stripped, 

high  Z  ions) ,  the  collisional  ionization  cross-sections  are 

modified  by  only  a  few  percent.  For  low  Z  ions  with  a  large 


number  of  bound  electrons  (such  as  Neon-like  ion  with  10 
electrons  per  ion)  and  at  higher  plasma  density,  the 
modifications  to  cross-sections  and  rates  may  be  substantial. 
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TABLE  1 

Li-like  A  l  XI 

Electronic 

Trans i t i on 

:  1  s  23d  ,  3 D  -* 

1 s  24s  ,  2 S 

Energy  of 

Partial 

Scattering 

incident 

Outgoing 

Scattering 

Total 

Cross- 

► 

Electron 

Electron 

Angular 

Cross- 

Collision 

Section 

( Ryd ) 

( Ryd ) 

Momentum 

Sect i on 

Strength 

( nan  2 ) 

0 

0.21 24E-03 

1 

0 . 3074  E-02 

2 

0.8544E-02 

► 

78.17 

1.36 

3 

0.1811 E — 0 1 

0.7374 

1.13 

• 

4 

0.2543E-01 

E-01 

E-1  9 

f 

5 

0 . 7966  E-02 

6 

0.6874E-02 

7 

0 . 2589  E-02 

8 

0 .6871 E-03 

0 

0.1 335  E-03 

1 

0.2254E-02 

2 

0. 5444E-02 

3 

0. 1 232E-01 

116.57 

39.76 

4 

0. 2346E-01 

0.6738 

6 . 92 

t 

5 

0.4668E-02 

E-01 

E-20 

1 

6 

0.7656E-02 

7 

0.5692E-02 

8 

0 . 3073E-02 

9 

0.1431 E-02 

0 

0 . 8864  E-04 

1 

0 . 1  7  42E-02 

2 

0. 3829E-02 

3 

0 . 8646E-02 

0.6344 

4 . 90 

154.98 

78.17 

4 

0.211 6E-01 

E-01 

E-20 

5 

0.3589E-02 

6 

0.5929E-02 

7 

0 . 6469E-02 

8 

0.4824E-02 

9 

0. 3037E-02 

1  0 

0.1749E-02 

L. 
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TABLE  1  (Cont’d) 

Li-like  A l  XI 


Electronic  Transition:  1s23d, 


-*■  1 s  2  4s  ,  2S 


Energy 

of 

Partial 

Scattering 

Incident 

Outgoing 

Scattering 

Total 

Cross- 

Electron 

Electron 

Angular 

Cross- 

Collision 

Section 

( Ryd ) 

(Ryd) 

Momentum 

Section 

Strength 

( irao  2 ) 

0 

0.4249E-04 

1 

0.11 56E-02 

2 

0.2303E-02 

3 

0 . 4646E-02 

231  .79 

154.98 

4 

0.161 8E-01 

0.6011 

3.10 

5 

0.4792E-02 

E-01 

E-20 

6 

0.3121 E-02 

7 

0 . 4878E-02 

8 

0.5310E-02 

9 

0.4692E-02 

1  0 

0. 3645E-02 

1  1 

0.2622E-02 

0 

0.2228E-04 

1 

0.8410E-03 

2 

0 . 1 6l 2E-02 

3 

0.281 5E-02 

4 

0.11 79E-01 

0.5740 

2.23 

308.60 

231  .79 

5 

0.6131 E-02 

E-01 

E-20 

6 

0.2678E-02 

7 

0.3271 E-02 

8 

0.41 32  E-02 

9 

0.4474E-02 

10 

0 . 4225 E-02 

1  1 

0 . 3595E-02 

1  2 

0 . 2857  E-02 

1  3 

0.21 66E-02 

0 

0.1231 E-04 

1 

0.6489E-03 

385.41 

308.60 

2 

0.1231 E-02 

0.5635 

1  .75 

3 

0.191 2E-C2 

E-01 

E-20 

4 

0.8508E-02 

5 

0.6468E-02 

6 

0 .31 96  E-02 

7 

0.2688E-02 

8 

0.3081 E-02 

385.41 

308.60 

9 

0 . 3607E-02 

0.5635 

1  .75 

10 

0 . 3896  E-02 

E-01 

E-20 

1  1 

0. 3787E-02 

1  2 

0. 3373E-02 

1  3 

0.2787E-02 

1  4 

0.2269E-02 

15 

0.1 808E-02 
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TABLE  1  (Cont'd) 

Li-like  A l  XI 


Electronic 

Transition : 

1s23d,  3  D 

-»  1  s  2 4 s  ,  2S 

Energy  of 

Partial 

Scattering 

Incident 

Outgoing 

Scattering 

Total 

Cross- 

Electron 

Electron 

Angular 

Cross- 

Collision 

Section 

( Ryd ) 

(Ryd) 

Momentum 

Section 

Strength 

( Tran  2  ) 

0 

0.3989E-05 

1 

0.4329E-03 

2 

0.8266E-03 

3 

0.1151 E-02 

4 

0.4578E-02 

5 

0.5434E-02 

0.61 62 

1  .37 

539.04 

462.27 

6 

0.4044E-02 

E-01 

E-20 

7 

0.3049E-02 

8 

0.2501 E-02 

9 

0.2368E-02 

1  0 

0 . 2622E-02 

1  1 

0.2981 E-02 

1  2 

0.31  79E-02 

13 

0.3046E-02 

1  4 

0.2837E-02 

15 

0 . 2520  E-02 

0 

0.8388E-06 

1 

0.2789E-03 

2 

0.5441  E-03 

3 

0.7673E-03 

4 

0.211  8E-02 

5 

0. 3348E-02 

769.47 

692.66 

6 

0.3663E-02 

0.3276 

5.10 

7 

0.3565E-02 

E-01 

E  —  2 1 

8 

0.31  09E-02 

9 

0.2435E-02 

1  0 

0.1  930E-02 

1  1 

0.181  8E-C2 

1  2 

0.201 6E-02 

13 

0.221 3  E-02 

1  4 

0.2437E-02 

15 

0.251  5  E-02 
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TABLE  2 


Li-like  ML  XI 

Electronic  Transition:  1s22s,  2S  -*■  1s23S,  2S 


Energy- 

Incident 

Electron 

(Ryd) 

of 

Outgoing 

Electron 

(Ryd) 

Total 

Collision 

Strength 

Scattering 
Cross- 
Section 
( iran  2 ) 

252.58 

1  .36 

0 .9397E-01 

2.23E-19 

378.19 

126.97 

0.9070E-01 

1 . 44E-1 9 

503.80 

252.58 

0.8874E-01 

1 . 05E-1 9 

755.02 

503.80 

0.8633E-01 

6.84E-20 

1 006.25 

755.02 

0.8484E-01 

5.05E-20 

1257.47 

1006.25 

0.8380E-01 

3.99E-20 

1759.91 

1 508.69 

0.8243E-01 

2.80E-20 

251 3-57 

2262.35 

0.81 23E-01 

1 . 93E-20 
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TABLE  3 

Li-like  AX,  XI 

Electronic  Transition:  1s23s,  2S  -*•  1s23d,  2D 


Energy 

Incident 

Electron 

(Ryd) 

of 

Outgoing 

Electron 

(Ryd) 

Total 

Collision 

Strength 

Scattering 
Cross- 
Section 
( iTan  2 ) 

86.53 

1  .36 

0.4201 E  +  00 

2.91 E-1 8 

129.12 

43.95 

0. 4223E+00 

1  . 96E-1 8 

171.71 

86.53 

0.4232E+00 

1  .  4  7  E  —  1  8 

256.88 

171.71 

0. 4228E+00 

9 . 85  E-1 9 

342.06 

256.83 

0.4221 E+00 

7.38E-1 9 

427.23 

342.06 

0.4211 E+00 

5.90E-19 

597  .58 

512.40 

0 . 41 91 E+00 

4.20E-1 9 

853.10 

767.93 

0.41 65E+00 

2.92E-19 
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TABLE  4 
Li-like  A  £  XI 

Electronic  Transition:  1s23s,  2S  ■+  1s23d,2D 


Excitation 


Temperature 

(OK) 

Rate  Coefficient 
cm  !sec' 1 

3 . 50E+05 

1 . 39E-1 2 

5 . 80E+05 

1 . 82E-1 0 

8 . 1 0E+05 

4.35E-10 

1 . 00E+06 

5 . 97E-1 0 

1 . 20E+06 

6.78E-1 0 

3 • 50E+06 

7.30E-10 

5 . 80E+06 

7 • 32E-1 0 

8 . 1 0E+06 

6 . 35E-1 0 

1  . 00E+07 

5.63E-10 

3 . 50E+07 

4.80E-10 
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TABLE  5 


Li-like  A l  XI 

Electronic  Transition:  1s22s,2S  -*•  1s23s,2S 

Excitation 


Temperature 

(OK) 

Rate  Coefficient 
cm 3  sec" 1 

3 . 50E+05 

3.28E-20 

5 . 80E+05 

1 . 64E-1 3 

8 . 1 0E+05 

3.44E-12 

1 . 00E+06 

1 .21 E-l 1 

1 . 20E+06 

2.1 5E-1 1 

3.50E+06 

3.1 8E-1 1 

5 . 80E+06 

8.94E-1 1 

8 . 1 0E  +  06 

9.51 E-1 1 

1 . 00E+07 

9.18E-1  1 

3 • 50  E+07 

8.37E-1 1 
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TABLE  6 


Li-like  hi  XI 

Electronic  Transition:  ls23d,2D  -►  1s2iis,2S 


Temperature 

(OK) 

Excitation 
Rate  Coefficient 
cm  3sec~  1 

3 . 50E+05 

1 .07E-1 3 

5.80E+05 

7.93E-12 

8 . 1 0E+05 

1 .65E-1 1 

1 . 00E+06 

2.11 E— 1 1 

1  . 20E+06 

2.31E-1  1 

3  •  50E+06 

2.42E-1  1 

5 . 80E+06 

2.06E-1 1 

8 . 1 0E  +  06 

1 .62E-11 

1 . 00E+07 

1 • 33E-1 1 

3 • 50E+07 

1 .OHE-1 1 
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